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Abstract 

Studies of DD correlations for a large sample of events containing fully and par- 
tially reconstructed pairs of charmed D mesons recorded by the Fermilab photo- 
production experiment FOCUS (FNAL-E831) are presented. Correlations between 
D and D mesons are used to study heavy quark production dynamics. We present 
results for fully and partially reconstructed charm pairs and comparisons to a re- 
cent version of Pythia with default parameter settings. We also comment on the 
production of Y>(3770) in our data. 



1 Introduction 



Heavy quark production continues to present itself as a challenge to our under- 
standing of the strong interaction. While Quantum Chromodynamics (QCD) 
provides a theoretical framework for our understanding and perturbative QCD 
can be applied to some aspects of heavy quark production, other aspects re- 
main elusive and cannot be described without including a variety of non- 
perturbative effects. This is especially true for charm production, where per- 

* See http://www-focus.fnal.gov/authors.html for additional author information. 
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turbative QCD calculations involve large uncertainties and non-perturbative 
effects play a significant role in modeling physical observables. Until we achieve 
a fundamental understanding of the strong interaction, accurate models that 
are able to reproduce properties of the strong interaction — such as heavy quark 
production — are crucial for our understanding of this fundamental force. 

In this paper, we present new results from FOCUS (FNAL-E831) on charm- 
pair correlations between D and D mesons. Charm-pair correlations have re- 
ceived considerable theoretical attention [1,2,3,4,5,6], and have been studied 
in both hadropro duct ion [7,8,9,10,11,12] and photoproduction [13,14] experi- 
ments. We present our photoproduction results by comparing data distribu- 
tions to predictions from a recent version of a Monte Carlo based on the Lund 
Model [15], which includes non-perturbative effects that have been shown to 
be important in charm production. We select default settings for charm pho- 
toproduction in the Monte Carlo to facilitate comparisons with theoretical 
predictions and results from other experiments. 



2 Experimental method 

The data for our studies of DD correlations were recorded by the FOCUS 
experiment during the 1996-1997 fixed-target run at the Fermi National Ac- 
celerator Laboratory. The experiment ran with a photon beam 1 and a spec- 
trometer that was upgraded from a previous photoproduction experiment, 
E687 [16]. The FOCUS spectrometer had a target that consisted of four BeO 
target elements for most of the recorded data 2 . A vertex detector, which was 
located in the target region, had a total of 16 planes of silicon strip detectors. 
Four of the planes were interleaved with the BeO target elements, and 12 were 
located downstream of the target. Tracks that were reconstructed in the ver- 
tex detector were linked to particle tracks that were found in five multiwire 
proportional chambers. Particle momenta were determined by measuring the 
deflection of tracks in two analysis magnets of opposite polarity, and parti- 
cle identification was accomplished using measurements from three multicell 
threshold Cerenkov counters, details of which are described elsewhere [17]. 

Here we describe, for the first time, the candidate-driven algorithm that was 
used to collect a large sample of ~ 7000 pairs of fully reconstructed charmed 
mesons. The sample consists of pairs of D mesons: D + D~, D + D°, D°D~, 

1 The photon beam was produced from the bremsstrahlung of secondary electrons 
and positrons with an endpoint energy of ~ 300 GeV. The average photon energy 
for the recorded data was ~ 180 GeV with a width of ~ 50 GeV. 

2 Early in the run a few different targets were used, and less than 5% of the charm- 
pair data were recorded with Be (instead of BeO) target elements. 
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Fig. 1. (a) Normalized D invariant mass vs. normalized D invariant mass distribu- 
tion, and (b) a fit to the normalized D invariant mass after sideband subtraction 
(described in the text), (c) Invariant mass of the recoil D in the partially recon- 
structed charm-pair sample (the mass of charged D candidates is lowered by 3.74 
MeV/c 2 to match the D° mass distribution). The yield is a sum of individual yields 
for the three decay modes, (d) A 2 distributions for right-sign (filled triangles) and 
wrong-sign (open circles) combinations for partially reconstructed charm pairs. 

and D°D°. For this paper, we considered the decay modes D° — > K~tt + , 
D + — > K~it + it + , D° — > K~-k + 'k + 'k~ , and charged-conjugate modes. The 
algorithm considered all combinations of two, three, and four charged tracks 
to find a combination that could be associated with the decay of a single 
D meson, and a second combination of tracks that could be associated with 
a second D decay vertex in the same event. The successful reconstruction 
of two D vertex candidates was followed by the reconstruction of a primary 
interaction vertex, particle identification cuts, and detachment cuts for the D 
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vertices relative to the primary vertex. The goal was to achieve low background 
levels for each decay mode using a minimum number of cuts. 

The first step of the candidate-driven algorithm considers all pairs of two-, 
three-, and four-track combinations in an event. Each combination of tracks 
represents a possible D decay. For each track the algorithm considers all pos- 
sible combinations of charged K or 7r assignments such that the assignments 
are consistent with the decay of a charged or neutral D meson. A particular 
combination of tracks and the associated particle assignments is referred to as 
a D candidate. The mass of each D candidate is calculated using the measured 
track momenta, and is required to fall within a wide range of 1.6-2.4 GeV/c 2 . 
To select events with a D and a D, the kaons for the two D candidates are 
required to have opposite charge. 

The second step is vertex reconstruction. The goal is to find a pair of IP- 
decay vertices that can be associated with a primary interaction vertex, and 
to find all other tracks in the event that can be associated with that primary 
vertex. This part of the algorithm starts by performing a vertex fit for each 
D candidate. The tracks for each D candidate are required to form a vertex 
with confidence level greater than 1%. Pairs of D candidates that satisfy the 
confidence level cut are subjected to two additional vertex cuts. The first 
cut requires that the momentum vectors of the two candidates intersect with 
a confidence level greater than 1%. The second cut rejects background by 
rejecting pairs of D candidates for which the reconstructed daughter tracks 
for both D candidates form a single vertex with confidence level greater than 
0.1%. This rejects background events in which tracks for both candidates all 
come from a common vertex. The final phase of the vertex reconstruction 
treats the two D candidates as seed tracks to find the primary vertex. Vertex 
fits are performed by including the two seed tracks as well as combinations 
of all other tracks in the event. As many tracks as possible are added to the 
primary vertex as long as the confidence level is greater than 1%. 

Pairs of D candidates that survive the vertex reconstruction are subjected 
to particle-identification cuts, which are based on measurements from three 
multicell threshold Cerenkov counters. The Cerenkov algorithm [17] calculates 
four likelihoods that correspond to the four hypotheses (electron, pion, kaon, 
proton) that are considered for each charged track. The algorithm produces a 
X 2 -like variable Wi = — 2 ln(likelihood), where % is the index used to represent 
each hypothesis. For the kaon in each D candidate, we require that the kaon 
hypothesis is favored over the pion hypothesis by more than a factor of exp(0.5) 
by requiring W n — Wk > 1.0. For the pions in each D candidate we apply a 
pion consistency cut, which requires that no particle hypothesis is favored over 
the pion hypothesis with a AW = W n — W m i n greater than 5, where W min is 
the Wi with the smallest value. 
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Fig. 2. (a) Invariant DD mass for D + D~ and D°D mass combinations for back- 
ground-subtracted FOCUS data (open circles), Pythia 6.203 (solid line), FOCUS 
data with iV pr i mary = 2 cut (filled triangles), and Pythia 6.203 with iVp r i mary = 2 
cut (dashed line). The inset shows the invariant DD mass that we obtain after 
applying additional cuts, such as cuts that remove events with energy deposited 
in the electromagnetic calorimeters, (b) Number of tracks assigned to the primary 
vertex for background-subtracted FOCUS data (data points with error bars) and 
Pythia 6.203 (solid line) normalized to the number of DD pairs in data with 
-^primary > 2. (c) Acj) and (d) pf of the DD pair for background-subtracted FOCUS 
data with -/V pr i mary > 2 (open circles), E687 data (filled triangles with offset to show 
error bars) normalized to FOCUS data, and Pythia 5.6 (solid line). 

After applying particle-identification cuts, we impose cuts based on the signifi- 
cance of detachment (tjot) between each D candidate and the primary vertex. 
We calculate £/<j£ by using the measured value of £, the distance between the 
D decay vertex and the primary vertex, and dividing by the associated error 
og. The cuts for Ijai range from Ijag > 1 to if on > 4 depending on the 
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decay mode, whether the .D-decay vertex is located between target elements 
(for which background levels are low) or in target material, and whether a D 
candidate can be associated with a D* decay. 




Fig. 3. Correlations for fully reconstructed DD pairs with iVprimary > 2: (a) A<f>, 
(b) pi of the DD pair, (c) rapidity difference (yo — and (d) invariant DD 
mass for background-subtracted FOCUS data (data points with error bars), Pythia 
6.203 after detector simulation and data analysis cuts (solid line), and Pythia 
6.203 parent distributions without acceptance or resolution effects (dashed line with 
arbitrary normalization) . 

Figure la shows the DD signal that we obtain after all of the aforementioned 
cuts have been applied to the FOCUS data. Figure la shows the normalized 
D invariant mass 3 M n (D) opposite the normalized D invariant mass M n (D). 

3 The normalized mass, M n (D) = AM/ctm, is defined as the difference between the 
reconstructed mass and the central value of the D + or D° mass distribution divided 
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Figure lb shows a Gaussian fit to M n (D) over a linear background after apply- 
ing a background subtraction procedure that is used to determine the number 
of pairs of charmed D mesons in the FOCUS data. The procedure consists of 
performing a sideband subtraction and fit for one normalized mass distribution 
by selecting entries in the signal and sideband regions of the other normalized 
mass distribution. In Fig. lb we plot M n (D) by assigning unit weight to D 
candidates with a reconstructed mass in the signal region (± 2er about the 
central value of the D + or D° mass of the candidate), and a weight of — 1 / 2 to 
candidates with mass in the two 4-8<r sideband regions 4 . The DD yield that 
we obtain from our fit is 7064 ± 119 (statistical error). 

In addition to our study of correlations between pairs of fully reconstructed D 
mesons, we study correlations between two D mesons where one D is fully re- 
constructed and the other is kinematically tagged by a slow pion coming from 
the decay D* + — > ir + D°. In these decays, the D° need not be reconstructed, 
and therefore we refer to this sample of charmed D mesons as partially re- 
constructed charm pairs 5 . The reason for including this sample in our studies 
of correlations is that charm-pair correlations can be studied over a larger 
kinematic range compared to the fully reconstructed sample. 

For partially reconstructed charm pair events [14] we begin by considering all 
two-, three-, and four-track combinations for the fully reconstructed D (recoil 
D) in an event. We consider the decay modes D° — > K~n + , D + — > K~it + -n + , 
D° — > K~it + it + it~ , and charged-conjugate modes. A candidate-driven algo- 
rithm uses the recoil D candidates to find the primary vertex, requiring the 
vertex confidence level to be greater than 1%. The same Cerenkov particle 
identification criteria used for the fully reconstructed charm-pair sample (see 
above) are applied to the recoil D candidates. However, a more restrictive 
detachment cut of tjot > 5 is applied to all three decay modes in the partially 
reconstructed charm-pair sample. Figure lc shows the invariant mass distri- 
bution, which includes all three decay modes, with a total of 782 630 ± 1600 
candidates satisfying the selection criteria. 

The next step in the analysis treats each track that is assigned to the primary 
vertex (excluding the recoil D) as a slow-pion candidate from the decay D* + — > 
n + D°. The momentum of the track is multiplied by 13.8 to approximate the 
momentum of the D* + 6 . If the charge of the slow pion is the same as the 

by the reconstructed- mass error gm, which is calculated for each D candidate. 

4 An equivalent approach to determine the DD yield is a fit to M n (D) after selecting 
signal and sideband regions for M n (D). Using this approach we obtain a DD yield 
of 7126 ± 120 , which is consistent with the yield mentioned in the text. 

5 The partially reconstructed sample consists of D* + D~, D* + D°, D°D*~, and 
D + D*~ pairs. 

6 Due to the low Q value of the D* decay, the momentum of the soft pion ap- 
proximates the momentum of the D* when multiplied by the inverse of its energy 
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charge of the kaon from the recoil D, then the combination of the slow pion 
and recoil D is designated as a right-sign combination. Otherwise, it is a wrong- 
sign combination. This assignment of right- and wrong-sign combinations is 
used for background subtraction. 

A double subtraction method is used to reduce backgrounds. First, to handle 
non-charm background, a sideband subtraction is applied to recoil D candi- 
dates. A Gaussian fit is applied to the invariant mass distribution for each of 
the three decay channels. Entries in the 4-8(7 sideband regions are subtracted 
from those in the ± 2a peak region by using a weight factor of — 1 / 2 . Sec- 
ond, the assignment of right- and wrong-sign combinations is used to subtract 
wrong-sign background from right-sign combinations. To avoid distortion of 
the wrong-sign background we exclude all slow-pion candidates that can be 
associated with a D* decay involving the recoil D. This anti-D* cut is im- 
posed by excluding combinations of slow pions and recoil D-mesons that have 
a mass difference, m(D*) —m(D), in the range 0.142-0.149 GeV/c 2 . To further 
enhance the selection procedure, a maximum cut of 4 (GeV/c) 2 is applied to 
A 2 = (pi r) + 13.8 * pW) 2 + (pW + 13.8 * pW) 2 , where pt\p { y ] and p£\pW 
are transverse momentum components of the recoil D and slow pion, respec- 
tively. This cut enhances the selection of signal since genuine events balance 
A 2 (see reference [14] for more details). This is shown in Fig. Id, which shows 
a prominent excess of right-sign combinations close to A 2 = compared to 
the wrong-sign background. After applying the double subtraction and the A 2 
cut, we obtain a sample of 75 160 ± 1040 partially-reconstructed charm pairs. 



3 DD correlations 

For our study of correlations between pairs of fully reconstructed D mesons, 
we compare FOCUS data to predictions from a Monte Carlo based on the 
Lund Model. The Monte Carlo consists of a Pythia 6.203 [15] generator with 
default settings, and detector simulation algorithms for the FOCUS appa- 
ratus. The Monte Carlo generator produces charm events using a tree-level 
photon-gluon fusion process applied to beam photons and target nucleons. 
We use default options for charm photoproduction in the generator (instead 
of using a Monte Carlo tuned to match our data) to facilitate comparisons 
with theoretical predictions and results from other experiments. In this paper, 
we also compare our results to previously published charm photoproduction 
results from experiment E687 [14]. 

To improve comparisons between data and model predictions based on photon- 
gluon fusion, we eliminate our lowest multiplicity events by requiring a mini- 
fraction, which is « 13.8. 
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Fig. 4. Charm-pair correlations for the partially reconstructed charm-pair sample: 
(a) A(j), (b) pi of the DD pair, (c) rapidity difference (yp — Vp)-, an d (d) invariant 
DD mass for background-subtracted FOCUS data (data points with error bars) 
and Pythia 6.203 after detector simulation and data analysis cuts (solid line). The 
A<p distribution for fully-reconstructed charm pairs (open circles with error bars) is 
included in (a) after accounting for resolution broadening (see footnote 10). Pythia 
parent distributions (dashed lines with arbitrary normalization) are included in (c) 
and (d). The slight distortion (dip) at the peak of the Ay distribution in (c) is 
caused by the anti-D* cut described in the text. 

mum number of particles assigned to the primary interaction vertex. We define 
-^primary as the number of particles assigned to the primary vertex. With this 
definition, iVp rimary has a minimum value of two since it includes the D and D 
mesons (each charm meson counts as a single particle) in addition to charged 
tracks assigned to the primary vertex. To eliminate our lowest multiplicity 
events we require an A^ primary > 2 cut. The cut eliminates features observed 
in data that are not present in Pythia 6.203. This is illustrated in Fig. 2a, 



10 



which shows the background-subtracted 7 invariant DD mass for mass com- 
binations with a net charge of zero (D+D- and D°D ) for FOCUS data, 
and for Pythia events that have passed through a software simulation of the 
FOCUS detector and have survived the event selection procedure described 
earlier in this paper. The mass distribution has an enhancement near threshold 
that is not present in Pythia. This enhancement is evident for events with 
-^primary = 2, especially when we apply additional cuts that remove events 
with energy deposited in electromagnetic calorimeters (see inset in Fig. 2a). 
The enhancement seems to arise from the diffract ive production of -0(3770) 
decaying to DD, and will be the subject of a future paper (additional in- 
formation can be found in conference proceedings [18]). Another significant 
difference between data and Pythia is shown in Figs. 2a and 2b. This is the 
excess of N pr - unary = 2 events in data compared to Pythia, some of which 
can be attributed to the production of ■0(3770). By eliminating iVp rimary = 2 
events, we get fairly good agreement for the iV primary distribution in Fig. 2b, 
which shows the histogram for Pythia (solid line) normalized to the number 
of DD pairs in the data (data points with error bars) with A^ primary > 2. By 
eliminating the iV primary = 2 bin the agreement between data and Pythia is 
significantly improved (a slight excess of events with A^ primary = 3 persists in 
the data). 

Previous studies [7,8,9,10,11,12,13,14] of charm-pair correlations have pre- 
sented distributions for pf(DD), the transverse momentum squared of the 
DD pair, and A0, the azimuthal angle between the D and D momentum 
vectors in the plane transverse to the beam direction. These distributions are 
significant, since pf(DD) = and A0 = 7r radians in leading-order QCD, 
where the charm-quark pair is produced back-to-back. In QCD these distri- 
butions are broadened by NLO corrections and non-perturbative effects, as 
illustrated in references [3] and [4]. Photoproduction results from E687 [14] 
have been compared to results from NLO calculations [1] and Pythia version 
5.6 [19]. The E687 comparisons between data and Pythia 5.6 are reproduced 
in Figs. 2c and 2d, but with a different normalization to match FOCUS data 
(shown as open circles with error bars). The figures show good agreement be- 
tween FOCUS and E687 data, and a significant discrepancy between data and 
Pythia 5.6. 



7 The background subtraction procedure assigns unit weight to DD candidates in 
the signal region in Fig. la (± 2a about the center of the distribution), a weight 
of — 1 /2 to candidates in the single D and single D sidebands (four regions defined 
as ± 2a about the D axis and ± A-Sa about the D axis, and ± 2a about the D 
axis and ± 4-8<r about the D axis), and a weight of + 1 / \ to candidates in the 
four regions where both the D and D candidates are 4-8a away from the center 
of the distribution. The weight factor of + 1 / \ accounts for the over-subtraction of 
the single- Z) and single- D backgrounds and the subtraction of random combinatoric 
background. 
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Agreement between FOCUS data and the more recent Pythia 6.203 is sig- 
nificantly better, but minor discrepancies persist. Figure 3 shows comparisons 
for A<p, Pt(DD), rapidity difference defined as Ay = y D — yjy, and invariant 
DD mass, M(DD). FOCUS data are plotted as data points with error bars. 
Pythia parent distributions (dashed lines) are shown without acceptance or 
resolution effects, so that parent distributions can be compared to the distri- 
butions that are obtained for Monte Carlo events that have survived detector 
simulation, event selection and analysis cuts (solid histograms). 

Figure 3a shows good agreement for A0. There is an enhancement in the first 
A(f) bin, which is not present in Pythia and may suggest the presence of 
an additional production mechanism. There is good agreement for p^(DD) in 
Fig. 3b, except that the data tend to have slightly larger values of pf(DD). 
Compared to Pythia 5.6, the agreement between data and Pythia 6.203 
for A<fi and pf(DD) is significantly better. Some of the improvement can be 
attributed to a larger value for the intrinsic transverse momentum of the in- 
coming partons, referred to as the kr kick 8 , but a number of other Pythia 
modifications that affect these distributions have also occurred over time. Fig- 
ure 3c shows fairly good agreement for Ay 9 , but also shows significant ac- 
ceptance losses for |Ay| > 1 (acceptance losses are less severe in the partially 
reconstructed charm-pair sample). Acceptance losses are also significant for 
large values of M(DD) in Fig. 3d, but here there is a discrepancy between 
data and Pythia for smaller values of M(DD) where the acceptance is good. 

Figure 4 shows results for the partially reconstructed charm-pair sample, com- 
paring data (asterisks with error bars) to Pythia 6.203 (solid lines). For A<p 
(see Fig. 4a) we also include a comparison to the distribution that we obtain 
for fully reconstructed charm pairs after accounting for resolution broadening 
effects 10 . This shows that the two samples are in agreement, and that the en- 
hancement that we observe in the first A<p bin for fully reconstructed charm 
pairs (see Fig. 3a) disappears due to resolution broadening and selection cuts 
applied to the partially reconstructed charm-pair sample. Figures 4a and 4b 
are both affected by resolution broadening (the effects are reproduced by our 
Monte Carlo), and the agreement between data and Pythia 6.203 is good. 
As before, the data tend to have slightly larger values of p^(DD). 

8 A value of < k\ > = (1 GeV/c) 2 was introduced with Pythia version 6.135, 
while previous versions had a value of < k\ > = (0.44 GeV/c) 2 

9 The agreement between data and Pythia improves slightly for D mesons with 
larger values of £/crg, however a more restrictive i/cr^ cut also reduces the number 
of charm-pair events that are available for correlation studies. 

10 The A(j) distribution for the fully-reconstructed sample is obtained by taking the 
momentum vector of the D or D in an event and treating it as the momentum of a 
D* that decays isotropically to a D° and a pion. The pion momentum vector is then 
used to determine Acj) as is done in the analysis of partially reconstructed charm 
pair events. 
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In Figs. 4c and 4d we show results for Ay and M(DD), and include Pythia 
parent distributions (dashed lines) to show how acceptance losses in this sam- 
ple compare to acceptance losses in the fully-reconstructed sample (see Fig. 3). 
The partially reconstructed charm pairs are less affected by acceptance losses, 
and thus extend the kinematic range of our correlation studies. The Ay dis- 
tributions in Fig. 4c show good agreement, while the M(DD) distributions in 
Fig. 4d exhibit a mismatch between data and Pythia that is similar to the 
mismatch that is observed in Fig. 3d. 



4 Conclusions 

We have extracted two large samples of photoproduced charm-pair events for 
studies of correlations between D and D mesons. The first sample consists 
of more than 7000 fully reconstructed DD pairs. The second sample con- 
sists of over 75 000 partially reconstructed charm pairs, where one D meson 
is fully reconstructed and the other is tagged by a slow pion coming from a 
D* decay. For the fully reconstructed sample we impose an A^ primary > 2 cut 
to eliminate our lowest multiplicity events, while the partially reconstructed 
sample has an implicit cut of Af primary > 2 due to the presence of the slow 
pion. The significance of the N priinary cut is that it improves our comparisons 
to model predictions based on photon-gluon fusion by eliminating low multi- 
plicity events in which we observe the production of ?/>(3770) decaying to DD 
pairs. The -^(3770) events, which are not included in Pythia, appear to be 
produced diffract ively, and will be the subject of a future paper. 

The FOCUS results on charm-pair correlations presented in this paper are in 
good agreement with previous measurements from experiment E687, which 
displayed significant discrepancies compared to an older version of Pythia 
(version 5.6). Comparisons of FOCUS data to a more recent version of Pythia 
(version 6.203) are significantly better, due to changes in parameters that 
affect the modeling of photon-gluon fusion. One notable change that im- 
proves the agreement with data is that the intrinsic transverse momentum 
(k T ) of incoming partons was increased from < k\ > = (0.44 GeV/c) 2 to 
< k\ > — (1 GeV/c) 2 . Although minor discrepancies persist when FOCUS 
data are compared to Pythia, the modeling of heavy quark photoproduction 
is fairly good for correlations between D and D mesons. 
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